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INTRODUCTION 


Three Javelin rocket flights from Fort Churchill, Canada have provided 
detailed data on electron and proton precipitation in different types of 
auroral forms. These rockets named TWINS 1, TWINS 2 , and p/» were launched 
in 1967 and 1968. A thorough understanding of auroral phenomena, requires 
a high time resolution study of electron and proton energy spectra and 
pitch angle distributions with the ability to separate spatial and temporal 
effects. The WINS 1 and 2 data furnish this .necessary information and 
in addition study such things as the characteristics of the backscattered 
electrons and the transit times of the different energy auroral particles. 
Analysis of these data will give greater insight into the details of 
auroral production mechanisms. Under this contract these rocket data as 
well as pertinent satellite data measured in the magnetosphere will be 
studied to search for the mechanisms responsible for the auroral phenomena. 


ACCOMPLISHMENTS 

During the first year of the contract the TWINS data on precipitated 
and backscattered auroral electrons was compared with backscattered fluxes 
derived from a theoretical coulomb scattering analysis due to Walt, MacDonald, 
and Francis. It was found that the measured backscattered fluxes in the 
energy range of .5 keV to 2 keV often exceeded the theoretically predicted 
fluxes. This discrepancy can be explained by an electric field parallel 
to the magnetic field with a total potential drop of about 1 KV. Although 
this parallel field can substantially alter the characteristics of the pre- 
cipitating auroral electrons, it is not of sufficient magnitude to act as 
the original source of these particles. A reprint from the Journal of Geo - 
physical Research covering the details of this portion of the work is en- 
closed in this report. 
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During the second contract year we have pursued the details of the theo- 
retically predicted hackscattered electron fluxes. Our work has been directed 
toward combining the backscatter code of Walt et al. ( 1969 ) which is accurate 
down to energies of about 500 eV with a code devised by Banks and Nagy which 
was originally used to predict photoelectron fluxes both upgoing and down- 
going in and above the neutral atmosphere. Dr. Nagy and Dr. Banks have 
worked to extend their code up to energies of 500 eV where it has been 
meshed with the higher energy results of Walt et al . This combined pro- 
gram is now able to predict the entire hackscattered electron energy 
spectrum from auroral energies of 20 KeV down to a few eV. The combined 
computer code will also predict the height emission profiles of the im- 
portant auroral optical emissions and can be used eventually to derive 
the auroral ionosphere created under a given arc. 

A typical run of the combined program would be done as follows. 

Any given auroral spectrum covering the full energy range (a few eV to 
10* s of KeV) can be used as input. The portion of the spectrum greater 
than 500 eV is input into the Walt program. This program uses Coulomb 
scattering and the Bethe continuous energy loss formula to diffuse the 
particles in energy and pitch angle. ■ It calculates the new distribution 
of degraded primaries as a function of energy (above 500 eV), pitch angle, 
and altitude. The degraded primary fluxes are then used directly with 
the neutral density profile and the appropriate ionization cross sec- 
tions to calculate the secondary electrons produced at each altitude. 

These secondary electrons are then used as input for the low energy 
portion of the program. This part is computed using the Nagy-Banks 
portion with the initial energy spectrum below 500 eV and the secon- 
daries as a function of altitude as the input. The Nagy-Banks portion 
of the program then calculates the transport and energy degradation of 
these secondary electrons, and the resulting distribution is used to 
calculate the different emission and ionization profiles. 
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The details of the combined auroral electron code are discussed 
in the attached preprint which has been submitted to the Journal ' of 
Geophysical Research . This paper shows many examples of incident 
auroral spectra and "monoenergetic" electron distributions. The 
spectral degradation, backscatter spectra, ionization rates, and 
optical emission profiles for these spectra are also shown. This 
preprint gives an excellent overview of the scope of the new auroral 
code which represents the bulk of the second year's contract effort. 

During the fourth quarter of the second year the TWINS II data 
were examined for possible further comparison with the computational 
results of the newly developed program. ' In particular, a period 
toward the end of the flight in which the rocket penetrated the at- 
mosphere in a nearly vertical direction was used to correlate the 
actually measured degrading electron spectrum with that calculated 
by the code. This approach proved to be unfeasible, however, because 
the temporal variations in the incident auroral spectrum were suffi- 
ciently great that they obscured the variations due to atmospheric 
degradation. 


NEW TECHNOLOGY 

The instruments discussed here were developed under previous 
NASA grants. Consequently there has been no new technology developed 
by this contract. 
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CONCLUSION 


During the second year of this contract a significant new 
tool has been added to the study of auroral physics. The newly 
developed auroral electron code is now able to predict the scat- 
tering and energy loss of electrons covering the complete auroral 
energy spectrum from .5 eV up to several MeV. The code can be 
used very effectively in input-output experiments in which the 
incident auroral electron spectrum is measured by a satellite 
passing overhead and the resulting ionization or optical emissions 
are measured on the ground. The code will need to be continually 
refined to improve its accuracy and to e'xpand its capabilities, 
but even in its present form it enables the modeling of the inter- 
action of auroral electrons with the atmosphere to be done in a 
detailed way not previously possible, and consequently offers a 
major contribution to the study of the aurora and the ionosphere. 
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Energy spectra and pitch angle distributions of auroral electrons have been measured in a 
premidnight multiple arc auroral display by a Javelin rocket containing two identical payloads 
that separated in flight. The rocket (code-named Twins 2) was launched from Fort Churchill, 
Canada, at 0459 UT on March 2, 1968, and covered an altitude range up to 800 km. The elec- 
tron energy spectra between 40 ev and 20 kev show a 'continuum' spectrum with a superim- 
posed energetic peak. The center energy of the peak was-observed to shift from 10-12 kev 
over the arcs to 2-3 kev between the arcs. This spectral structure is shown to be similar to the 
inverted ‘V’ structure reported by other investigators. The in flight separation of the two pay- 
loads allowed investigation of spatial versus temporal effects in the auroral precipitation. In 
one interval of the flight an arc was observed to be moving northward with a velocity of 
0.6 km/sec. Calculated backscattered spectra are compared with those actually measured. Good 
agreement was observed for electron energies above about 1.5 kev, but below this energy the 
backscatter ratio was observed to be — 100% . Several explanations for this unusually high ratio 
are considered, and a likely possibility is shown to be a parallel downward-directed electric 
field between ~150 and 500 km. 


Measurements of auroral particle fluxes, their 
spectra, and pitch angle distributions have been 
the object of numerous satellite and rocket in- 
vestigations. A recent review of pertinent mea- 
surements has been given by Pavlikas [1971], 
The general features of auroral particle pre- 
cipitation are now well known, although the 
details of the origins of the particles and of 
their interaction with the atmosphere are far 
from understood. Several investigators have re- 
ported ‘anomalous’ results in measurements of 
auroral particles. For example, Mozer and 
Bruston [1966a, 6] reported energetic protons 
traveling upward from mirror points well within 
the atmosphere and anticorrelations between 
electron and proton fluxes. A recent study by 
Reme and Bosqued [1971] reported observa- 
tions of anticorrelated proton and electron fluxes 
and proton pitch angle distributions peaked 
along field lines. These observations were inter- 
preted as indicating the presence of electric 
fields parallel to magnetic field lines in the 
upper ionosphere with magnitudes ranging up 
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to a few hundred millivolts/meter. Furthermore, 
Mozer and Bruston [1967] reported direct ex- 
perimental evidence by means of rocket-borne 
electric field probes of such parallel electric 
fields in the upper ionosphere. 

Choy and Arnoldy [1971] have observed an 
anomalously high reflection coefficient or back- 
scatter ratio for auroral electrons with energies 
below 1 kev. These results are indeed surprising 
in view of the relatively low altitude of the 
measurements (<150 km, R. L. Arnoldy, pri- 
vate communication, 1972). 

Chappell [1968] has compared calculated 
and measured auroral electron backscatter ratios 
using data from the Twins 1 twin payload 
sounding rocket [ Westerlund , 1969]. The theo- 
retical backscatter fluxes were calculated by 
means of a Fokker-Planck treatment of atmo- 
spheric Coulomb scattering. Discrepancies be- 
tween calculated and measured backscatter 
were suggestive of atmospheric heating, a 
parallel electric field below the rocket, or a 
combination of both effects. 

In this paper we report observations of pre- 
cipitating (30° < a < 60°) and backscattered 
(120° < a < 150°) auroral electron fluxes in 
the energy range 40 ev to 20 kev over a series 
of IBC-2 auroral arcs. The measurements were 
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made with the Twins 2 payloads, a pair of 
identical payloads that were launched on an 
Argo D-4 Javelin sounding rocket and were 
separated in flight. Five distinct regions of 
particle precipitation are identified in terms 
of total energy input, and differing particle 
spectra characteristics are shown for each re- 
gion. Variations in spectra near the edge of an 
arc are shown to resemble the ‘inverted V’ spec- 
tral structure [ Frank and Ackerson, 1971], 
Finally, measured and calculated backscatter 
spectra are compared, and an attempt is made 
to explain the observed discrepancies. 

Description of the Experiment 

A detailed description of the Twins 2 sound- 
ing rocket experiment and the auroral condi- 
tions at launch may be found in O’Brien and 
Reasoner. [1971]. Briefly, the payloads were 
launched from Fort Churchill, Canada, on 
March 2, 1968, at 0459 UT into an IBC-2 
northward-moving multiple arc structure. The 
ground magnetometers recorded a 100-y nega- 
tive bay at the time of the flight. 

The rocket carried two identical payloads, 
Twins 2A and 2B, which were separated in 
flight with a relative velocity along the flight 
direction at separation of approximately 4 
m/sec. The relative and absolute velocities per- 
pendicular to . the magnetic field lines were 
1 ± 0.2 m/sec and 585 ± 6 m/sec, respectively. 
The flight azimuth was 37° east of true north. 
Auroral electron spectra were measured over an 
altitucle range of 200-800 km, and the payloads 
covered L values from L ss 9 to L ss 13. 

Each payload carried identical detector com- 
plements, consisting of two wide range electron 
spectrometers, three fixed-energy differential 
electron detectors, and a large-area Geiger tube. 
Of relevance to this study were the two electron 
spectrometers, code-named Specs [ O’Brien et 
al., 1967]. The Specs detectors were electrostatic 
deflection spectrometers employing channel elec- 
tron multipliers and were sensitive to electrons 
with 40 ev < E < 20 kev. The Specs were 
laboratory calibrated both by an electron beam 
and by placing a calibrated ro Ni source over the 
apertures, thereby establishing the relative effi- 
ciencies of the upward- and downward-looking 
instruments. The upward-looking Specs scanned 
a pitch angle range of 30°-60° in ~5 9 incre- 
ments, and the downward-looking Specs simul- 


taneously scanned a range of 120°-150°. The 
Specs required 3.6 sec to complete a measure- 
ment cycle of ah electron energies. Payload 
aspect information was provided by a flux gate 
magnetometer. 

Experimental Results 

An 1 overview of the auroral conditions en- 
countered during the Twins 2 flight is best 
shown by a plot of the total energy deposited, 
obtained by integrating over the measured 
spectra. Figures 1 and 2 show the values of 
the directional electron energy deposited and 
backscattered as calculated from the measure- 
ments of Twins 2A and 2B, respectively, for 
electron energies >400 ev and pitch angles of 
30° < a < 60° and 120° < a < 150°. The 
payloads encountered two distinct periods of 
enhanced precipitation at 240-350 sec and again 
at 670-710 sec. We can associate these periods 
of enhanced precipitation with the visible au- 
roral arcs: The precipitated fluxes were iso- 
tropic to within a factor of 2-3 over the pitch 
angle range 30° < a < 60°. When it is assumed 
that the isotropy extended over the upper hemi- 
sphere, the total energy deposited was therefore 
~20 ergs/cm 2 sec. This agrees wefl with the 
photometrically measured brightness of the au- 
roral forms at launch of 10-20 kR. 

The gross features of Figures 1 and 2 are in 
good agreement, as would be expected, since the 
maximum payload separation was 700 meters. 
The backscattered flux measured by Twins 2A 


ALTITUDE (KM) 423 744 769 480 

L 9.5 10.7 11.7 12.5 



Fig. 1. Total directional incident and back- 
scattered electron energy flux measured by Twins 
2A. The measurements were computed by in- 
tegrating the measured electron spectra, appro- 
priately weighted, over the energy range 400 ev 
to 20 kev. 
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ALTITUDE (KM) 423 744 769 480 

L 9.3 10.7 11.7 12.5 



Fig. 2. Total directional incident and back- 
scattered electron energy flux measured by Twins 
2B. The measurements were computed by in- 
tegrating the measured electron spectra, appro- 
priately weighted, over the energy range 400 ev 
to 20 kev. 

is apparently less than that measured by Twins 
2B, but this was due to a failure of two of the 
six channels of the downward-looking Specs on 
Twins 2A. Therefore, part of the electron spec- 
trum was not measured by this detector. The 
general characteristics of both incident and 
backscattered particle spectrums as measured 
by both payloads were similar throughout the 
flight, thus increasing our confidence in the data 
presented herein. We have elected to concen- 
trate on spectrums measured by Twins 2B 
because of the partial failure of the Twins 2A 
downward-looking Specs. 

A closer examination of Figures 1 and 2 
reveals that the decrease in precipitation at 
T = 350 sec was observed 6 ± 3 sec earlier 
by Twins 2B (the following payload) than by 
Twins 2A (the leading payload). This means 
that the arc observed from 250 to 350 sec was 
moving northward more rapidly than the pay- 
loads and effectively overtook and ran ahead 
of them. Considering the payload absolute hori- 
zontal velocity of 0.585 km/sec and the separa- 
tion at this time of 217 meters, the arc is cal- 
culated to be moving northward with a velocity 
along the payload flight direction (37° east of 
north) of 0.61 km/sec. Flux enhancements later 
in the flight (e.g., T = 630 sec) do not display 
any time delays between detection by the two 
payloads, although the Specs detectors with 
their 3.6-sec cycle time would not have been 
able to resolve the time delay of ~2 sec ex- 
pected if the payloads encountered a stationary 


or slowly moving spatial structure. It is there- 
fore possible that the payloads initially observed 
a northward-moving arc, which subsequently 
slowed and was again encountered by the north- 
ward-moving payloads. The unique capability of 
the Twins payloads to distinguish temporal and 
spatial effects demonstrated above was also in- 
strumental in the discovery of short-duration 
field-aligned electron bursts reported previously 
by O’Brien and Reasoner [1971], 

Figures 3-8 show precipitated and backscat- 
tered electron spectra measured at various times 
in the flight by Twins 2B. These spectra repre- 
sent regions of differing precipitation intensity 
at T — 281, 414, 479, 630, and 670 sec. These 
times are indicated , by the arrows in Figure 2. 
Each spectrum represents an average over 3 
Specs cycles, or about 10 sec. The periods were 
selected to insure that there were no fluctua- 
tions in the data, which exceeded the usual 
variations due to statistics. Each point is an 
average of 147 samples of 10 msec each with 
the 147 points divided (by the Specs switching 
sequence) into three groups of forty-nine 10- 
msec samples. 

Figure 3 shows the incident, or precipitated, 
spectrum at T — 281 sec. The fluxes in units 
of el/cm 2 sec ster ev are plotted as a function 
of electron energy in electron volts. This spec- 
trum was measured at an altitude of 585 km 
and represents pitch angles of 32°-56°. The loss 
cone (referenced to the 100-km level) at 585 



ENERGY (eV ) 


Fig. 3. Incident electron energy spectrum mea- 
sured at T = 281 sec. For reference with regard 
to the time of measurement for this figure and 
Figures 4-9, see Figure 2. 
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Fig. 4. Backscattered electron energy spectrum Fig. 6. Incident and backscattered electron en- 
measured at T — 281 sec. . ergy spectra measured at T — 479 sec. 


km is 64° , and hence all of the observed incident 
electrons would theoretically mirror below 100 
km and thus be subject to severe atmospheric 
scattering and loss. All downgoing electrons 
measured throughout the flight were in the loss 
cone; hence no. information on the trapped- 
particle population was acquired. The group 
of data points at each energy represents the 
spread in flux at a given energy for different 
pitch angles betwen 32° arid 56° and indicates 
the variation of energy spectrum with pitch 
angle. This variation is a factor of 2 to 3 maxi- 
mum. The error bars for the points are also 
indicated. The errors are dominated by uncer- 


tainties in the instrument calibration of ±20%. 
Errors due to counting statistics are negligible, 
typically <5%. The precipitating electrons 
shown in Figure 3 have a peak in the energy 
range of 6-17 kev, which contains the bulk of 
the total precipitating’ energy. 

Figure 4 shows the electron spectrum that 
was backscattered from the atmosphere below 
the payload as measured simultaneously at T — 
281 sec by the downward-viewing Specs. The 
data format is the same as that for Figure 3. 
The pitch angle range measured was 124°-148°. 
Note that the backscattered spectrum exhibits 
a small peak in the range 5^15 kev, which is 
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Fig. 8. Incident and backscattered electron en- 
ergy spectra measured at T = 670 sec. 

almost an order of magnitude lower in intensity 
than the precipitated flux. , 

As the payloads passed from the region of 
intense precipitation into one of intermediate 
precipitation ( T — 414 sec), the spectra 
changed to those shown in Figure 5. Here we 
have shown only the sketched incident and 
backscattered spectra . without including the 
data points showing the spread in pitch angle. 
This is done here and in subsequent presenta- 
tions for the sake of brevity and to afford 
easier comparison of spectra. 

The incident spectrum at T = 414 (solid 
line) shows two peaks, a prominent peak in the 
range 2-9 kev ancTa less prominent one in the 
range 500 ev to 1.2 kev. This last peak, al- 
though it was not observed consistently through- 
out the flight, is nonetheless a real feature 
representing enhancements relative to the ‘con- 
tinuum’ in two adjacent energy channels with 
a difference of greater than 3 standard devia- 
tions. Note that, as the payloads passed out of 
the region of enhanced precipitation (Figure 3), 
the peak shifted down from 6-17 kev to 2-9 
kev. The backscattered spectrum (dashed, line) 
shows no evidence of a peak. Note also that the 
flux of backscattered electrons with E < 600 
ev is comparable to or greater than the corre- 
sponding flux of precipitating electrons. We will 
comment further on this surprising backscatter 
ratio in a later section. 

Figure 6 shows the incident and backscattered 
spectra in the region of low-intensity precipita- 


tion at T = 479 sec. The . peak in the incident 
spectrum has shifted to the energy range 1.5-4 
kev and has diminished in differential flux level 
by a factor of 2 from the previous spectrum. 
This precipitated spectrum measured between 
the intense precipitation regions is similar to 
the background continuum discussed by Wester- 
lund [1969] using data from Twins 1, a similar 
rocket payload fired i year earlier. Here, the 
backscattered spectrum shows no peak in 
energy, unlike the backscattered spectrum in 
the intense precipitation region (f = 281 sec, 
Figure 4). Again; the backscatter ratio is 
—'100% for energies less than 500 ev. 

Following the low-intensity region the rocket 
encountered another region of intermediate in- 
tensity precipitation with, the incident electron 
energy spectrum showing a characteristic up- 
ward shift in the energy of the peak to the 3- 
to 12-kev range. Figure 7 shows representative 
spectra from this region measured at T = 630 
sec. The backscattered spectrum has a slight 
peak in this same energy range. In this case the 
backscatter ratio is —100% for energies less 
than 1.2 kev. 

Finally, the payloads again encountered a 
region of intense precipitation (Figure 8) at 
T = 670 sec. The energy of the peak is in the 
3- to 16-kev range, and the incident spectrum 
is similar to that observed during the first en- 
hanced period at T = 281 sec. This adds sup- 
port to the thesis advanced earlier that the two 
regions of intense precipitation were one and 
the same arc. The backscatter peak is again 
prominent, and the backscatter ratio is —100% 
for energies below 1 kev. 

From this series of spectra it is clear that 
changes in the total energy deposited were 
caused by shifts in the energy of an energetic 
peak that appeared, to be superimposed on a 
much more stable continuum spectrum. During 
this Twins 2 flight, unlike the case in the Twins 
1 flight [ Westeriund , 1969], the peak never 
completely disappeared, although it was con- 
siderably diminished in both intensity and aver- 
age energy in the period of T — 460 to T = 
500 sec. 

A higher time resolution study of this peak 
shift effect is shown in Figure 9. Here three 
successive incident spectra at T = 351, T = 
366, and T = 405 sec are plotted to display 
the spectra] change over the time period of a 
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Fig. 9. A series of three incident electron 
energy spectra measured in the period T = 351 
to T = 405 sec showing the shift in the energy 
of the peak as the payload exits the region 
over an auroral arc. 

sharp decrease in total precipitating intensity 
(Figure 1). The shift of the energy of the peak 
to progressively lower values is clearly evident. 
We will discuss this further in a later section. 

Perhaps the most surprising feature of these 
spectra is the very high, — -100%, backscatter 
ratio for electron energies less than about 1 kev. 
The expected ratio computed from Coulomb 
scattering is usually of the order of 10%. In 
the following section we consider the back- 
scatter phenomena in detail, presenting calcu- 
lations of expected backscattered fluxes for com- 
parison with those actually observed. 

Comparison of Experimental and 
Theoretical Backscatter Results 

To examine the interaction of the precipitated 
auroral electrons with the atmosphere, we have 
used the measured incident electron spectrum 
as input for a theoretical computer code that 
calculates the expected backscattered electron 
spectrum. This calculated backscattered spec- 
trum can then be directly compared with the 
backscattered spectrum measured simultane- 
ously by the downward-viewing Specs detector 
on the Twins payload. Any differences between 
the calculated and measured backscattered spec- 
tra can be interpreted in terms of the physical 
processes occurring in the interaction region. 

The theoretical calculations were made using 
the Fokker-Planck diffusion equation solution of 
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Walt et al. [1968]. This treatment includes 
angular scattering, continuous energy loss, and 
the effect of a converging magnetic field. The 
solution is in the form of a distribution function 
for the auroral electrons as a function of energy, 
pitch angle, and altitude. The calculations were 
made using differential cross sections given by 
the first' Bom approximation [ Mott and 
Massey, 1965] and the continuous energy loss 
formula of Bethe [1933]. These calculations are 
reasonably accurate (±10%) for isotropic pitch 
angle distributions and for electron energies 
down to about 500 ev. The computations done 
in this work are carried down to 1 kev. 

Figure 10 shows the measured incident (solid 
line) and backscattered (circles) spectra at T = ' 
281 in the first region of enhanced precipita- 
tion. The incident spectrum, which was assumed 
to be isotropic in pitch angle over the upper 
hemisphere, was used as input to the backscat- 
ter code. The dotted area represents the calcu- 
lated backscattered flux in the range of pitch 
angles of 126°-149° corresponding to the pitch 
angles measured by the downward-viewing 
Specs. The agreement between the calculated 
backscattered flux and the measured flux (cir-. 
cles) is reasonably good above about 1.5 kev. 
However, below 1.5 kev the measured back- 
scattered fluxes are roughly equal in magnitude 
to the incident fluxes ‘ (backscatter ratio of 


— INCIDENT -OBSERVED PA. 32'-56" 

• BACKSCATTERED -OBSERVED P.A. 124*-148" 

EJ BACKSCATTERED -CAECULATED PA 126*-149". «>||-0 



Fig. 10. Comparison of measured and calcu- 
lated backscattered electron spectra at T = 281 
sec, a region of intense precipitation. 
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- ~100%). The theoretical calculations, which 
employ Coulomb scattering only, predict a 
backseatter ratio more on the order of 10-20% , 
not 100%. In other words, the electrons of 
energy below 1-2 kev in this enhanced pre- 
cipitation region are either kept out of the at- 
mosphere preventing any loss or the incident 
primary particles are replaced by secondaries 
produced in atmospheric scattering. 

One obvious explanation for the 100% back- 
scatter below 1-2 kev is the presence of a 
parallel downward-directed electric field above 
the atmosphere with a total potential drop of 
about 1 kv in magnitude. This field would 
prevent electrons of less than 1 kev from enter- 
ing the atmosphere and would reflect them un- 
diminished in flux value. Other potential ex- 
planations also exist for this unusually high 
backseatter ratio. The production of multiple 
secondary electrons in the atmosphere by the 
energetic precipitating primaries is a possibility. 
Also, a severe modulation of the magnetic field 
by the auroral electrojet currents could possibly 
affect the particle pitch angles and cause reflec- 
tion. In anticipation of the results of a later 
discussion on the relative merits of these possi- 
bilities, we will examine theoretically the effect 
of a parallel downward-directed electric field of 
potential drop of 1 kv situated below the rocket 
and above the effective scattering layer of the 
neutral atmosphere. In this T = 281 case this 
field would be between 585 km and about 150 
km. 

For this T = 281 case we have again used 
the measured incident spectrum shown in Fig- 
ure 10 but have included a vertical potential 
drop $n ’of I kv just below the rocket. The 
effect of this parallel field is to alter the pitch 
angle distribution of the incident’ particles in 
such a way as to reflect particles whose energy 
component parallel to the field is less than 1 
kev and to reduce the energy of particles that 
traverse the field region by 1 kev. The Coulomb 
scattering, energy loss, and mirroring of these 
penetrating electrons are then calculated by the 
computer code, and the resulting backscattered 
spectrum is then transformed back through the 
potential drop with the corresponding change 
in pitch angle distribution and increase in 
energy. The results of this calculation are shown 
by the shaded area in Figure 10 for energies 
down to 1 kev. Note that the effect of the field 


is to increase the flux backscattered at the more 
nearly perpendicular pitch angles for electrons 
between 1 and ~10 kev in energy and also to 
cause the .100% backseatter below energies of 
a few kev. Qualitatively then, the calculations 
show that, with the addition of a parallel down- 
ward-directed electric field above the atmo- 
sphere to the general effects of Coulomb scat- 
tering, the predicted backscattered flux is in 
good agreement with the measured backscat- 
tered flux for all energies. 

Figure 11 shows the measured and calculated 
backscattered fluxes for the time period of very 
low intensity precipitation, T = 479 sec. As 
before, the measured incident spectrum is shown 
by the solid line and is used as the input for 
the calculation. The shaded area shows the pre- 
dicted range of backscattered flux for the pitch 
angles measured by the downward-viewing 
Specs. The measured backscattered electrons 
are again shown by the circles. In this calcula- 
tion we have not included a parallel electric 
field. One can see that the agreement above 1 
kev for Coulomb scattering is quite good. Recall 
that this set of spectra was measured between 
the auroral arc regions. Although no parallel 
electric field is required to obtain agreement 
between calculated and observed backscattered 
spectra above 1 kev in this low-intensity pre- 
cipitation region, one can see that below about 
500 ev the backscattered flux again approxi- 
mately equals the incident flux, indicating the 



Fig. II. Comparison of measured and calcu- 
lated backscattered electron spectra at T = 479 
sec, the time of the minimum precipitation. 
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possible need for a lower-magnitude parallel 
field with potential drop of about 500 volts to 
give agreement. 

Figure 12 shows the set of spectra represent- 
ing the second payload encounter with the in- 
tense precipitation region at T = 670 sec. 
Again, the shaded area represents the predicted 
backscattered fluxes with no parallel electric 
field, and the circles represent the measured 
backscattered fluxes. As in the first case (T = 
281), the agreement at higher energies (greater 
than about 1.5 kev) is good. However, as in the 
previous encounter with the intense precipita- 
tion region, the backscattered fluxes below about 
1.5 kev are comparable to the incident fluxes in 
contrast with the predictions of the Coulomb 
scattering calculations, which show the back- 
scattered flux decreasing with decreasing energy 
below 1.5 kev. As before, one can see that the 
addition of a parallel electric field of about l'-kv 
total potential drop below the rocket would 
cause an increase in the predicted backscattered 
fluxes at lower energies (less than about 6 or 8 
kev) with a backscattered flux equal to the 
incident flux below about 1.5 kev, in agreement 
with the measured fluxes. 

We have found then that in regions of aurora) 
precipitation the measured backscattered elec- 
tron fluxes agree quite well with the back- 
scattered fluxes predicted from Coulomb scat- 
tering alone for the higher energies greater 
than about 1.5 kev. However, for the lower 
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Fig. 12. Comparison of measured and calcu- 
lated backscattered electron spectra at T = 670 
sec, a period of intense precipitation. 
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energies, agreement between predicted and mea- 
sured backscattered fluxes for Coulomb scatter- 
ing alone is not sufficient but requires an effect 
such as a parallel downward-directed electric 
field between the rocket and the atmosphere 
with a total potential drop of about 1 kv in 
regions of intense precipitation and about 500 
volts in regions of less intense precipitation. 

Summary and Discussion 

The features of the measured precipitated 
and backscattered auroral electron spectra can 
be summarized as follows. 

1. The spectra contain two main compo- 
nents, a continuum spectrum that remains 
fairly stable throughout the flight and a peak 
in the 6- to 17-kev range, which is present over 
the auroral arcs but which shifts to lower 
energy and intensity as the payloads move 
away from the arc. This peak, when it is pres- 
ent, carries the bulk of the total energy de- 
posited. 

2. The incident electron fluxes at all energies 
appear to be isotropic in pitch angle to within 
a factor of 2 or 3 over the observed pitch angle 
range of 30°-60°. 

3. The backscattered fluxes for energies be- 
low about 1 kev appear to be isotropic within 
a factor of 2 or 3 and equal in intensity to the 
precipitated fluxes. Above about 1 kev the 
fluxes show a larger variation with pitch angle 
as expected from the effects of atmospheric 
scattering and the magnitude of the back- 
scatter ratio (j{E) up /J(E) down) at a given 
energy decreases with increasing energy. 

4. A theoretical backscatter calculation em- 
ploying Coulomb scattering, continuous energy 
loss, and magnetic mirroring and using the mea- 
sured incident spectrum as input is able to 
reproduce the measured backscattered spectrum 
for energies greater than 1.5 kev in intense 
regions of precipitation and for energies of 
greater than 500 ev in regions of less intense 
precipitation. The measured backscattered spec- 
tra below about 1 .5 kev, however, cannot be 
explained on the basis of Coulomb scattering 
alone. 

In the first feature above we mentioned the 
shift in energy of the peak in the incident elec- 
tron spectrum as the payloads move in and out 
of the auroral arc region. If this shift in peak 
energy were displayed on an energy-time spec- 
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trogram, it would resemble an inverted V struc- 
ture similar to that discussed by Frank and 
Ackerson [1971], The widths of the structures 
observed here (50-100 km) are in reasonable 
agreement with the ‘typical’ widths of 200 km 
discussed by Frank and Ackerson [1971], Exact 
agreement on these widths is not expected owing 
to uncertainties resulting from high velocities 
of the auroral forms relative to the rocket pay- 
• load velocities. 

The most interesting feature of the electron 
spectra is the 100% backscatter ratio for ener- 
gies less than 1.5 kev. As we mentioned in the 
previous section, there are at least three pos- 
sible explanations for this large ratio : secondary 
electrons that are produced by the high-energy 
primary electrons and that are scattered up- 
ward out of the atmosphere, a large variation 
in the ambient magnetic field caused by the 
nearby auroral electrojet that changes the 
electron pitch angles, or a downward-directed 
parallel electric field with a potential drop of 1 
kv above the atmosphere, which reflects lower- 
energy electrons before they undergo scattering. 

We discount the possibility that secondary 
electrons are responsible for the upward fluxes 
at energies between 40 ev and 1 kev. If these 
particles were secondaries, we would expect that 
the measured upward flux of these particles 
would be correlated with the measured incident 
primary energy flux of the particles greater 
than 400 ev. In Figure 13 we show the incident 
energy flux from Figure 1 along with the 40-ev 
precipitated and backscattered number flux. 
There is a good correlation between the 40-ev 
downgoing and upcoming flux, but there is no 
significant correlation between the precipitated 
primary energy flux and the 40-ev upcoming or 
backscattered flux. Thus the experimental evi- 
dence argues against secondary electrons as a 
source of these low-energy electrons. 

Theoretical evidence also argues against sec- 
ondaries as a source of the high-intensity low- 
energy influx. Most of the secondaries that are 
produced are below 300 km, and the calcula- 
tions of Nagy and Banks [1970] show that these 
electrons cannot escape the atmosphere because 
of the very high degree of scattering. However, 
these same calculations show that a 50-ev elec- 
tron produced at an altitude greater than 300 
km can escape significant scattering and travel 
freely upward. Cross sections for secondary 
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Fig. 13. A comparison of the total incident 
energy flux (from Figure 1) with the measured 
40-ev incident and backscattered flux. 


electron production by energetic primary elec- 
trons incident on N 2 have been measured by 
Opal et al. [1971], Taking the integrated N 2 
density above 300 km for a medium density 
atmosphere as 2 x 10 15 atoms/cm 2 [ Johnson , 
1965] and the precipitated energy as 3 ergs/cm 2 
sec ster with an average energy of 10 kev, we 
compute the expected upward secondary flux at 
40 ev from above 300 km of 0.2 X 10 5 el/cm 2 
sec ster ev. This is just about equal to the 
change seen in the 40-ev upward flux in the 
time period of T — 340-440 sec as the primary 
energy flux drops from about 3 ergs/cm 2 sec ster 
to 0.4 erg/cm 2 sec ster. Therefore, although 
secondary electron production above 300 km 
appears to contribute to the upflux of low- 
energy electrons, its contribution is only a minor 
part of the total observed upflux (0.2 x 10 5 
el/cm 2 sec ster ev out of about 1.5 X 10 E el/cm 2 
sec ster ev). 

Recent calculations by A. F. Nagy (private 
communication, 1972) have shown that the 
backscatter ratio below 40 ev can exceed 100% 
owing to backscattered and secondary electrons. 
However, above 40 ev the ratio falls off smoothly 
to the 20% value of 500 ev shown by the pres- 
ent calculations. 

As in the case of the secondary electrons, we 
can also rule out the possibility of magnetic 
field perturbations as a cause of the large back- 
scatter ratios. First, invariant calculations show 
that an increase of a factor of more than 3 in 
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the ambient magnetic field is required to cause 
a particle with a pitch angle of 30° at 585 km 
to mirror above the neutral atmosphere. Mag- 
netic field changes of this magnitude are not 
observed in the vicinity of aurorae. 

After discounting secondary electrons and 
magnetic field perturbations, the obvious ex- 
planation remaining is the presence of a paral- 
lel downward-directed electric field below the 
rocket of potential drop about 1 kv. A field of 
this sort would cause the incident auroral elec- 
trons with energies less than 1-2 kev to mirror 
above the region of significant scattering in the 
atmosphere. The backscatter calculations of the 
previous section have shown that the addition 
of such a field would give good agreement be- 
tween calculated and observed backscattered 
electrons. 

There are, however, certain crucial limitations 
on this electric field that must be met in this 
case. The first is that a total potential drop of 
~1 kv must exist between the rocket altitude 
and . the top of the scattering region (about 150 
km for 1-kev electrons). The second is that there 
appears to be no large fraction of this potential 
drop above 500 km, for there was not a system- 
atic large variation with altitude of the lowest- 
energy flux in the altitude range 500-800 km. 
Thus for the purpose of discussion we postulate 
a parallel electric field with a total potential 
drop of about 1 kev located between 150- and 
500-km altitude. 

Recent experimental and theoretical investi- 
gations have provided a possible source for the 
postulated electric field. Zmuda et al. [1966], 
Cloutier et al. [1970], and Choy et al. [1971] 
have inferred the presence of field-aligned cur- 
rents composed of very low energy electrons 
flowing out of the auroral ionosphere. Kindel 
and Kennel [1971] have theoretically investi- 
gated the effects of these currents and have 
concluded that ion acoustic and ion cyclotron 
instabilities, which are driven by the field- 
aligned currents, can occur under the conditions 
found in the auroral ionosphere. The role of 
the instabilities is to furnish an impedance to 
the flow of field-aligned currents and hence to 
cause a potential drop to develop. 

There is, however, a certain difficulty with 
the application of this mechanism to our par- 
ticular case. For the range of field-aligned cur- 
rent values expected just from flux balance in' 


the aurora studied here ( 1 0 s — 1 0 s el/cm 2 sec), an 
examination of Figure 7 of Kindel and Kennel 
[1971] shows that the first instability expected 
would be the 0 + of H + ion cyclotron mode and 
then only at altitudes greater than 1000 km. We 
require that the electric field, and hence the 
instability region; exist below 500 km. A pos- 
sible explanation for this apparent inconsistency 
is as follows. The critical current for instability 
is not only a function of altitude but also of 
the electron to ion temperature ratio. As the 
electrons are heated in relation to the ions, ion 
Landau and cyclotron damping become less im- 
portant in stabilizing the waves. An examina- 
tion of Figure 2 of Kindel and Kennel [1971] 
shows that, if the ratio of electron drift velocity 
to electron thermal velocity is initially above a 
certain value (VJV tb ~ 0.03 for TJ T t = 1), 
increasing the TJ T, ratio will cause waves to 
grow: Therefore, if the electrons in the 150- to 
500-km region could be heated sufficiently, the 
instabilities and hence the electric fields would 
exist and be confined to this region. Possible 
mechanisms for accomplishing this heating in- 
clude an energy transfer from secondary elec- 
trons by means of electrostatic waves (the 
‘bump on tail’ instability) or an upward trans- 
port of hot electrons from below. 

The above argument is necessarily qualita- 
tive, for data on crucial parameters such as 
temperature and drift velocity of the electrons 
contributing to the field-aligned currents are 
lacking. It is offered here as a plausible expla- 
nation for the required altitude structure of 'the 
postulated parallel electric field. Whatever the 
precise cause, the existence of the parallel down- 
ward-directed electric field appears necessary to 
adequately explain the data shown above: 

Conclusions 

The physical processes that describe the in- 
teraction of auroral electrons with the atmo- 
sphere appear to be more complex than just 
the Coulomb scattering of the incident primary 
electrons with a subsequent loss of energy. The 
comparison of the measured backscattered elec- 
tron spectra with spectra predicted using a 
theoretical scattering calculation has led to a 
discrepancy : for energies below about 1-2 kev. 
It was found that the very high ratio (~100%) 
of backscattered to incident fluxes for these 
energies couid most reasonably be explained by 
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a parallel downward-directed electric field that 
prevents these low r er-energy electrons from en- 
tering the atmospheric scattering region. This 
parallel field with potential drop of about 1 kv 
is thought to have its origin in wave-particle 
interactions^ in the turbulent auroral ionosphere. 
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ABSTRACT 


A new computational model has been created to describe the interac- 
tion of auroral electrons with the atmosphere. For electrons of energy 
greater than 500 eV continuous energy losses and small angle deflections 
are combined in a Fokker-Planck diffusion equation which computes energy 
spectrums over all pitch angles throughout the atmosphere. These fluxes 
are then used to determine the rates of secondary electron and degraded 
(E < 500 eV) primary electron production at all heights.. This informa- 
tion is used to compute upward and downward hemispherical fluxes in the 
energy range 0 to 500 eV taking into account discrete energy losses, 
large angle scattering and particle transport along magnetic field lines. 
The model has been used to compute energy spectrums, ionization rates, 
backscatter ratios and optical emissions associated with different inci- 
dent electron spectrums. For monoenergetic electrons of energy 2 keV 
and above the results obtained agree well with the work of Rees ( 19 69) 
and Rees and Maeda ( 1973 )* At lower energies the effects of trans- 
port and elastic collisions become progressively more important and 

the present results differ significantly both from the Rees and Maeda 
results and those obtained from the ideas of energy degradation. Finally, 
spectrums typical of the nighttime auroral oval and daytime polar cusp 
are used to obtain the altitude dependent fluxes, ionization rates and 
optical emissions. 
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I. Introduction 


There have been many previous attempts to provide an accurate des- 
cription of the way auroral electrons interact with the earth's upper 
atmosphere (see Rees, 1969 for a comprehensive review). The models re- 
sulting from these studies have become increasingly sophisticated in the 
past decade as new information about the various ionization and excita- 
tion processes has become available.. In addition, there has been a 
parallel improvement in our knowledge of the different aspects of the 
aurora itself, not only in terms of the nature of the incident electron 
energy and angular spectrums but also the different effects within the 
atmosphere; i.e. , the time and space variations of electron density and 
ion composition, the altitude variation of various optical emissions, 
the structure of the electron flux at different heights, and so forth. 

To a large extent many of the different models of auroral elec- 
tron behavior developed in the past have been based upon the ideas of 
energy deposition and, as a consequence, cannot be relied upon to accu- 
rately describe secondary ionization, dissociation, and other similar 
collisional interactions which depend upon the electron flux at energies 
less than 500 to 1000 eV. In an effort to improve this situation a new 
computational model has been developed to describe 'the penetration of 
auroral electrons into the upper atmosphere. This model incorporates 
both continuous and discrete energy losses while taking into account the 
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motions of electrons of all energies along the magnetic field, the var- 
ious effects of elastic and inelastic collisions, and the presence of 
secondary electrons arising from collisional ionization. 

The first widely used model of auroral electron penetration, de- 
veloped by Rees (1963), was based upon laboratory data for the range and 
energy deposition rate of electrons in air obtained by Griin ( 19 57 ) - Using 
an empirical expression based on the laboratory results, Rees computed 
the primary fluxes and the atmospheric ionization rate as a function of 
altitude for incident electron energies greater than U00 eV and for arbi- 
trary pitch angles. In a later study, Stolarski (1968) also used the Griin 
data to calculate the primary electron fluxes, secondary ionization rates 
and optical emissions. Monte Carlo techniques were first used by Maeda 
(1965) and Stadsnes and Maehlum (1965) for calculating the behavior of 
auroral electrons by considering continuous electron energy losses and 
angular scattering of many test electrons over a range of incident pitch 
angles. Certain results of the later Monte Carlo computations of Berger 
et al. (1970) were used by Rees (l 969) to create a new computational 
model but, for various reasons, the results of this hybrid method have 
been criticized (Shemansky, et al. , 1972) and further work has been 
necessary to overcome these difficulties (Rees and Maeda, 1973)* 

A different approach to the problem of auroral electrons was taken 
by Walt et al. , (1969) w ^° used a Fokker-Planck diffusion equation to 
describe simultaneously electron energy degradation, electron angular 
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scattering through collisions with the atmospheric gases and the effect 
of the earth's converging magnetic field upon the electron pitch angle 
distribution. In the method of Walt, et al. , the incident auroral elec- 
tron flux is described in terms of the pitch angle and energy distribu- 
tions. As the electrons pass into the atmosphere, analytic expressions 
for the energy loss (based on the Bethe formula) and angular scattering 
(based on the Born approximation) are used to compute the energy distri- 
bution and compute angular distributions (over both the upper and lower 
hemispheres) of the particles as a function of altitude. 

Although there are certain advantages associated with the work of 
Rees, Walt, et al. , Rees and Mae da, Berger and Seltzer and others, their 
methods share the problem that of any description of electron energy 
loss or angular scattering based upon continuous energy loss and small 
angle scattering becomes increasingly poor below 500 eV. Thus, while 
the behavior of electrons having energies greater than 500 eV can be 
accurately done using the expressions developed for a small fractional 
energy loss per collision, as one goes to energies below 500 eV the dis- 
crete nature of collisions in terms of energy loss and angular scatter- 
ing becomes of crucial importance. Although attempts have been made to 
compute separately the low (0 to 100 eV) energy portion of the secondary 
electron spectrum (Rees, et al. , 19 69 ), such an approach ignores the con 
tribution of the energy degraded primary flux and leaves unresolved ques 
tions about the details of the electron flux between 100 and 500 eV. In 
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addition such methods neglect the possibility of low energy electron 
transport along the magnetic field lines. 

Several years ago a new model of low energy electron behavior was 
introduced (Banks and Nagy, 1970; Nagy and Banks, 1970, 197l) in a study 
of the processes affecting the energy and space distributions of elec- 
trons created through atmospheric photoionization caused by absorption 
of solar EUV radiation. In this model two-stream continuity equations 
were used to describe electron fluxes moving into and out of the atmos- 
phere along magnetic field lines. Discrete energy losses and angular 
scattering were explicitly included in terms of appropriate energy de- 
pendent cross-sections. Although this, computational model was origi- 
nally intended to describe photoelectrons, it can be used equally well 
to describe the lower energy portion of the auroral electron energy 
spectrum. 

In the following sections a new model of auroral' electron behavior 
based upon continuous energy losses at high energies (E > 500 eV) and dis- 
crete .losses at low energies (E < 500 eV) is described. For high energies 
the Fokker-Planck diffusion approach of Walt et al. , (1969) has been ad- 
justed with a revised Born approximation (Chappell, 1968) to cover energy 
loss and angular scattering for electron energies ap low as 500 eV. Be- 
low 500 eV a discrete energy loss description based upon the method of 
Banks and Nagy (l970) and Nagy and Banks (1970) is used which takes into 
account the incident low energy particles, the energy degraded primary 
particles and the secondary electrons created through impact ionization. 
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The computer code resulting from the joining of the two methods gives 
complete energy spectrums of auroral electrons throughout the atmosphere 
in terms of the upgoing and downgoing electron fluxes. From these elec- 
tron fluxes it is possible to compute ionization rates for the different 
neutral gases, the rates of various optical excitations, and other perti- 
nent quantities. 

II. The Computational Model < 

Discrete excitation and ionization processes involving outer shell 
electrons dominate electron energy losses at low energies where the cross 
sections, while not known in complete detail, can be fairly well repre- 
sented through simple analytic expressions described later. Above about 
500 eV, however, progressively greater uncertainty is present with regard 
to specific processes involving both inner and outer shell electrons of 
the atmospheric gases. Further, as the electron energy increases, the 
fractional incident electron energy loss per collision decreases and the 
discreteness of the various losses becomes less important in describing 
the passage of electrons through the atmosphere. For the present model 
the joining point between the discrete and the continuous descriptions 
of electron energy loss has been taken as 500 eV. While somewhat arbi- 
trary, this choice reflects a compromise between computer data storage 
capacity and the accuracy of the continuous loss description at electron 
energies lower than 1 keV. (A more complete discussion of this point is 
given later.) In the following paragraphs descriptions are given of the 
continuous and discrete models. 
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Continuous Losses: Electron Energies Greater than 500 eV : 


To describe the behavior of electrons having energies greater than 
500 eV the Fokker-Planck diffusion equation approach of Walt et al. , 
(1969) has been used. In this model a magnetic field geometry with 
B oc R-3 is assumed to account for the effects of magnetic mirroring. The 
elastic scattering cross sections with the atmospheric gases are obtained 
by calculating the deflection of the incident electron by the Coulomb 
field of a charged scattering center. The atomic structure is relatively 
unimportant in this approximation and the differential scattering cross 
section can be written as 


o(B,8) - f lz (1 

4 m 2 c 4 3^ sin «/2 

where 6 is the scattering angle, Z is the scattering center charge, m and 
e are the electron mass and charge, 3 .= v/c, and c is the speed of light. 
The principal source of elastic scattering at high energies (E > 10 keV) 
arises from interaction with the nucleus; scattering collisions with the 
orbital electrons are relatively unimportant. At lower energies (l) must 
be modified to include the effects of orbital electrons. Such modifica- 
tion has been described by Chappell (1968) giving 


a(p,e) 


4m 2 c 4 


(z-F(e )) 2 


(1-e 2 ) 

3 4 


[l-B 2 sin 2 e/ 2 ] 
sin 4 9/2 


( 2 ) 


where F(fi) is an atomic scattering factor which introduces the shielding 
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effect of the orbital electrons on the scattering of lower energy 

(E < 10 keV) primaries. Values for F(9) are tabulated in the International 

Tables for X-ray Crystallography (1962). 

For the energy loss of the incident electrons, scattering by orbital 
electrons dominates and the Bethe formula is used in the form 


dT 

dx 


S njin (T me 2 J 2.7( T+2) /2I.) 
m 2 c 4 g 2 j J J 


(3) 


where T is the electron kinetic energy in rest mass units, I. is the aver- 

• J 

age excitation energy of the j-th species target atom (see Table l) and 
nj is the number density of bound electrons in the j-th atmospheric gas. 

Using equations (2) and (3), Walt, et al. , (1969) developed a 
Fokker-Planck diffusion equation to describe the passage of an incident 
electron flux through the atmosphere. The computations involve numerical 
solutions for the steady state electron distribution function with cal- 
culations being made of the electron flux as a function of altitude, 
energy and pitch angle. The input parameters for the model include the 
incident electron energy spectrum (E > 500 eV) and the initial pitch 
angle distribution. 

% 

As auroral electrons move into the atmosphere they are continually 
degraded in energy and eventually must either be backscattered out of the 
atmosphere or, pass through the computational threshold of 500 eV into the 
lower energy regime where the Bethe formula approximation of continuous 


8 



energy loss cannot be maintained. The energy degraded electrons are pro- 
duced at all altitudes in both upgoing and downgoing directions with a 
characteristic energy spectrum extending from 250 eV up to 500 eV with 
the precise energy dependence of the production rate at a given altitude 
reflecting the energy loss processes acting at energies above 500 eV. 
Although the altitude variation of the degraded primary electron pro- 
duction rate is proportional to the energy dissipation rate of the 
primary flux, it is necessary to note that these electrons are free to 
travel along the magnetic field to deposit their energy and create opti- 
cal excitations in regions separated from their point of origin. 


The incident high energy electron flux is also responsible for 
collisional ionization and the creation of secondary electrons. Like the 
degraded primary electrons, the secondary electrons are free to travel 
along field lines undergoing various inelastic and elastic collision pro- 
cesses. To model the production rate of the secondary electrons created 
through ionization impact collisions an analytical expression obtained 
from the work of Opal, et al., (l97l) has been used. If S(E ,E ) is the 

p b 

partial cross section describing the creation of a secondary electron in 
the range energy E g , E g + dE s by an electron originally possessing kinetic 
energy E^, the results of Opal, et al., can be approximated by the expres- 
sion 


S(Ep,£ s ) 


JL gj( E p) 

tS 1+ (E s /E) 2 


oo 
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where o^(Ep) is the ionization cross section at energy and E is an 
empirically obtained energy degradation coefficient. Values of E for 
N 2 , O 2 and 0 are given in Table 2. 


To model the ionization cross section as a function of electron 
energy the series of analytical expressions developed by A.E.S. Green and 
his collaborators have been used (Stolarski, personal communication, 1972. 
See also Green & Stolarski, 1972) • From this work the total electron 
impact ionization cross section can be written as 


Ci(E ) - 


Qn A 




V 


dW 


( 5 ) 


where q Q = 6. 51 x lO -1 ^ cm 2 eV 2 , I is the ionization potential of a par- 
ticular state of the final ion, A is an empirical parameter, W is the en- 
ergy loss of the primary particle, e = Ep/W and fi, y and v are chosen to 
permit ( 5 ) to fit experimental ionization cross-sections. 

The integration limits of ( 5 ) are chosen to accommodate the defini- 
tions of degraded primary and secondary electrons. In an ionization colli 
sion two electrons share what remains of the original incident electron 
kinetic energy (Ep-l) between them. The electron having kinetic energy 
in the range (E -l)/2 to (E -i) is referred to as the degraded primary 

lr }r 

electron while the electron having kinetic energy in the range 0 to 
(E -l)/2 is the secondary electron. If the primary electron loses an 
energy W, it follows that the secondary electron has. an energy W-I. A 
sketch of the way secondary and primary electrons are distributed in 
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energy as a result of a flux of electrons of energy E p is shown in 
Figure 1. 

Closed expressions corresponding to (5) can be obtained when v is an 
integer (see Table 3 for parameter values for N 2 , C> 2 and 0) : 

0: v = 1 Oi(E p ) = G(A-B) 

N 2 : v =2 ai(Ep) = G(A-2B+C) 

0 2 = v = 3 Oi(E p ) = G(A-3B+3C-D) 

where 
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Using the data of Table 3> ionization cross sections for the atmos- 
pheric gases can be computed from threshold to large, but non-relativist ic, 
energies. 

From the foregoing comments it follows that in addition to the inci- 
dent low energy (E < 500 eV) flux at high altitudes, there are two sources 
for the low energy flux within the atmosphere which arise from the passage 
of high energy (E > 500 eV) electrons through the neutral gases; i.e., the 
production of secondaries and degraded primaries. The altitude and energy 
dependent production rates for secondaries and degraded primaries can be 
obtained from (4) once the spectrum of the E > 500 eV electrons is known 
and a model atmosphere has been adopted. To illustrate this, the E < 500 eV 
production rate obtained for a typical nighttime auroral spectrum (discussed 
in a later section) is shown in Figure 2 as a function of electron energy 
for several altitudes. Owing to the form of the incident electron spectrum 
the simple symmetry shown in Figure 1 is absent since the production of 
secondaries dominates the production of degraded secondaries. 

Discrete Losses: Electron Energies Less Than 500 eV: 

For energies less than 500 eV the fractional energy loss per colli- 
sion becomes large and the effects of large angle elastic scattering 
arising from incident electron-orbital electron interactions cannot be ig- 
nored. Electron flux calculations for the low energy portion of the spec- 
trum are made using the two stream continuity equations described by Banks 
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and Nagy (1970), Nagy and Banks (1970) and Stolarski (1972). In this 
method the transport of electrons along the magnetic field lines is 
described in terms of the hemispherical fluxes of two electron streams $ + 
(upwards) and §" (downwards). These fluxes are defined as 

$ + (e,z) = j*0( E,z,n)an (6) 

Uppe r 
Hemisphere 

*"(E,z) = j 0 (E,z,n)dfJ (7) 

Lower 

Hemisphere 

where E,z,n) is the electron flux per unit energy and solid angle at 
the altitude z. 

The continuity equations for $> + and.$ _ can be derived by consider- 
ing: (a) inelastic collisions (energy degradation) with forward and 

backscattering and (b) elastic scattering (no significant energy degrada- 
tion) with forward and backscattering and (c) the production of electrons 
in a given energy range. When it is assumed that gravity, electric field 
parallel to magnetic field lines, and the effects of magnetic field con- 
vergence can be neglected, the steady state continuity equations for $ + 
and can be written as 


d$ + 

ds 


-1 

(cos 8) 


2 n k [ 


Os* + Pe k Oe k 


K 


(cos 0) Zj n k p e 


k k * - 


+ q/2(cos 9) + q + /(cos 0) 


( 8 ) 
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k k ,+ 
n kPe o e $ 


d$ ' 


-1 


ds (cos 0) 


1 <^STi) X 


+ q/2(cos 9) + q"/(cos 0) 


( 9 ) 


where 

_ k _ v k 
°a “ °ai 

q + (E,s) = E n k (s) E j (E')a^(E' -♦ E)$’(E' ; s) 

e>E. 

■ + (1 - p£i<E')a£j(E'- E) # + (E',s)| 


q-(E,s) = E D k (s) E |pE J (E-)a^j(E'-.E) t + (E',s) 
k J 

€>E 

+ (1 - Paj( E ')<^j(E'-»E)$-(E%s)J 

and 


s, 

$ (E,s) , 
$“(E,s), 
n R (s), 

Pe'W- 

O e h E ) r 


distance along a magnetic field line (positive outward), 
electron flux upward along s 
electron flux downward along s 
kth neutral species number density 

electron backscatter probability for elastic collisions 
with the kth neutral species 

electron total scattering cross section for elastic 
collisions with the kth neutral species 
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q(E,s), electron production rate in the range E to E + dE 
due to ionization processes; ie , the production of 
secondary electrons, the energy degraded primaries 
falling in the range 250 to 500 eV and, if necessary, 
the presence of photoelectrons created through the 
absorption of ultraviolet radiation. 


q 4 , electron production in the range E to E + dE due to 
cascading from higher energy electrons undergoing 
inelastic collisions. 


p ., electron backscatter probability for collisions with 

the kth neutral species resulting in the jth inelastic 
process . 



inelastic cross section for the jth excitation of the 
kth neutral species. 


To take into account pitch angle effects an average pitch angle 
(cos 9) has been introduced into (8) and ( 9 ); a pitch angle of 68° 
has been used in the calculations presented later. 


The cross sections and parameters used in the present calculations 
have been taken from Green and Stolarski (l972)and Stolarski (personal 
communication, 1972). The formulas describing ionization processes have 
been already given by equation (b) and ( 5 ) , and Tables 3 and 4 list the 
ionization cross section parameters used. 


The analytic form of the cross sections describing electron impact 
excitation of atomic and molecular states is taken from Green and 
Stolarski ( 1972) as 



j 


( 10 ) 
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where q Q = 6. 51 x 10“-^ cm^ eV^, W is the excitation energy and E is the 
incident electron energy. Table 5 lists the adopted values of the excita- 
tion parameters. 

The cross sections for elastic collisions for energies less than 
100 eV have been discussed previously by Nagy and Banks (l9T0). Table 6 
lists the values adopted here for the range 100 < E < 500 eV. The elas- 
tic cross sections for 0 were taken from Sunshine, et al. (1967) while 
the values of Fisk (1936) were used for Og and N^. For electron-thermal 
electron collisions the revised energy transfer rates of Schunk and Hays 
(l97l) have been adopted. . 

The elastic and inelastic backscatter probabilities are important 
parameters but, unfortunately, information about them is very sparse 
(see McDaniel, 1964; Mott and Massey, 1965; Kieffer, 1967; Lassettre, 

1969; Shyn, et al. , 1972, and Opal et al. , 197l) • The elastic back- 
scatter probabilities measured for Ng by Shyn, et al., (1972) were 
adopted for both Ng and Og elastic and excitation collisions. The back- 
scattering probabilities for 0 calculated by Fink and Yates (1970) have 
been used for both the elastic and excitation processes. For ionization 
processes it is assumed that the primaries continue in the same direction 
while the secondary electrons are distributed isotropically. 

To a large extent the success of the present method depends upon 
the matching of energy loss rates in passing from the continuous to the 
discrete regimes. A comparison of these energy loss rates is shown 
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for N 2 in Figure 3 * The dotted curve is derived from the discrete processes 
given in Tables 3 and 5 while the solid curve follows from the Bethe 
formula of equation (3)* At progressively lower energies the Bethe formula 
is first too large, then much smaller than the discrete losses. At the 
500 eV boundary used here the Bethe energy loss formula value is 2 . 9 % larger 
than the sum of the. discrete losses. For C >2 and 0 the discrepancies are 
l-bfi and 1.8$>, respectively. 

Method of Solution 

The computed auroral electron flux is obtained from sequential solu- 
tions to the Fokker-Planck diffusion equation and the low energy, two stream 
continuity equations. After the incident auroral energy spectrum and 
pitch angle distribution has been adopted, the electron flux for E > 500 
is computed as a function of pitch angle, energy and altitude. These 
fluxes are then used to determine, at each altitude, the rate of production 
of E < 500 eV degraded primary and secondary electrons. 

In the second step, the two stream continuity equations are solved 
using an implicit integration technique. The lower boundary condition for 
these equations, taken deep in the atmosphere, is § + (E,z-»0) = $~(E,z-+o) = 0 
for all energies. At the upper boundary $“(E,z = a) is specified by the 
incident low energy auroral spectrum. 

The low energy computations are started at '500 eV and work downwards 
taking into account energy cascading and scattering. Thus, although one 
can identify an initial degraded primary and secondary electron production 
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rate arising from ionizations caused by E > 500 eV electrons, further 
ionizations by E < 500 eV electrons are cumulatively computed such that 
the final spectrum cannot be segregated into a hierarchy of primary, ter- 
tiary, etc. electrons. The result obtained is the complete auroral energy 
spectrum which includes the effects of electron transport and multiple 
collisions. 

Once the auroral electron fluxes have been found, optical emissions, 
ionization rates and ionospheric heating rates are obtained using the 
different excitation and ionization cross sections described in the pre- 
ceding paragraphs. 

Ill Results 

As a first application of the present model we consider the behavior 
of nearly monoenergetic electrons having an isotropic pitch angle distri- 
bution in the downwards hemisphere. The results for these cases can be 
compared directly with the work of Rees (1969) and Rees and Maeda ( 1973) - 
Following the discussion of monoenergetic electron effects, results are 
given for more realistic electron fluxes typical of the nighttime auroral 
oval and daytime polar cusp. In all examples a vertical magnetic field 
(i = 90°) i s assumed, and a 1000°K thermospheric model of Banks and 
Kockarts (1973) has been adopted. 

Monoenergetic Electrons 

We assume that the incident electrons are isotropic in pitch angle 
over the lower hemisphere at an initial height of 985 km. The energy 
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distribution is Gaussian and has the form 

<i) = A 0 -(E-E o ) 2 /a, 2 (11) 

where the value a = 0.1 E Q is adopted for E Q = 0.42, 0.8, 2.0, 5-0, and 
10 keV. 

The normalization factor. A, for ( ll) has been chosen to agree with 
the unit incident flux normalization adopted by Rees (1963, 1969). If 
cp is the pitch angle dependent incident flux, the incident or, more com- 
monly, the perpendicular flux is defined through an integration over 
solid angle with the factor v cos 8; i.e. 

$ unit = (isotropic) (12) 

The fluxes plotted in the following figures represent hemispherical fluxes 
described by Equations (6) and (7) and the normalization demanded by (l2) 
has been chosen only for the purposes of comparison with other work. 

As energetic electrons penetrate the atmosphere low energy electrons 
in the range 0 to 500 eV are created through impact-ionization and excita- 
tion. The way in which these electrons are created has been briefly dis- 
cussed in Section I. Figure b shows the production rates of 0 to 500 eV 
electrons at several altitudes for an incident E Q = 10 keV flux. Although 
both secondaries and degraded primaries contribute to the total production 
rate over this energy range, cascading hides the simple distribution of 
secondaries and degraded primaries sketched in Figure 1. 
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Computations of electron flux as a function of altitude have been 
made for monoenergetic spectrums characterized by values of E^ in the 
range 0.42 to 10 keV. Figure 5 gives the vertical variation of $ + and $“ 
at several energies for E q = 10 keV. From these results it can be seen that 
although there is no incident [$“(z = 585 km)] flux present at the lower en- 
ergies, an upward (§ + ) flux is created through elastic and inelastic scat- 
tering. As shown, this backscattered, energy degraded, flux represents an 
outward flow of low energy electrons from the atmosphere. 

Another feature apparent in Figure 5 is the attenuation of the 
incident electron flux at the center energy E Q . At high altitudes these 
changes arise from magnetic mirroring while deeper in the atmosphere both 
pitch angle diffusion and energy loss are present. 

Information about the energy distribution of the downcoming electron 
fluxes at different altitudes is contained in Figure 6. As expected from 
the previous work of Rees, et al.(l969), Rees (1969), Rees and Maeda (l973)> 
Berger and Seltzer (1970) and others, the portion of the spectrum in the 
range 0 to 2 00 eV is steep, reflecting directly the way secondary electrons 
are produced (see Equation 4). This decline, coupled with the eventual 
rise in the flux when the presence of degraded primary electrons becomes 
important, has been discussed previously be Shemansky, et al. (1972) and 
Rees and Maeda ( 1973) - Although earlier calculations (see Rees, 1969, 
for example) could not properly bridge this transition region, the present 
results and those of Rees and Maeda (1973) give consistent agreement for 
energies at 2 keV and above. 
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A second interesting feature of Figure 6 is the tendency for the 
low energy portion of the electron spectrum to be relatively insensitive 
to changes in altitude. Such a behavior can be predicted in a crude way 
from theory. If we neglect transport and cascading, the rate at which 
secondary electrons are produced is given by no^(Ep) tf^Ep) where n is 
the neutral gas density, a^Ep) is the appropriate differential ioniza- 
tion cross section given by (k) and 0 ^(Ep) is the incident ionizing flux. 
Likewise, the loss of secondary electrons at the energy E g is given by 
n o^(E g ) 0 (E g ) where cr^(E g ) is the cross section for electron loss 
at the energy E g . If we equate production and loss we obtain the approx- 
imate formula 

*<*.> - £$$*!< v (13) 

which, since the neutral gas densities cancel, illustrates the altitude 
independence of the low energy flux. Of course, one cannot expect a 
simple relation like ( 13 ) to be strictly true; it only indicates the gen- 
eral trend of the results. However, we should also note that ( 13) pre- 
dicts the way that 0 (E g ) will depend upon the energy of the incoming 
flux. Since o^(Ep) is approximately proportional to E"-*- In E p for 
Ep ^ 1 keV, incident fluxes of progressively lower characteristic energy 
will produce proportionately greater numbers of low energy particles. 

With regard to upward electron fluxes, their energy spectrums are 
similar to the downgoing fluxes with a more pronounced minimum. Examples 
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of the upward flux spectrums are shown in Figure 7 forE Q = 10 keV. In 
the present models there is a substantial backscattering of the incident 
electrons arising from collisions and magnetic mirroring. A plot of 
downward/upward flux ratios is given in Figure 8 for e Q = 10 keV to illus- 
trate the altitude variation of the flux ratio and the net albedo of- the 
atmosphere to incoming electrons. The results for the low energy elec- 
trons show that although there are no incident low energy electrons 
(0 to 500 eV) , scattering within the atmosphere acts to create an upwards 
flux of substantial magnitude. At low altitudes the ratio 0 - /0 + is greater 
than unity, indicating the influence of the incident high energy particles. 
Above 300 km, however, 0 + is greater than 0" as large fluxes of low energy 
electrons leave the atmosphere. 

Plots of the ionization rates per unit incident flux for the mono- 
energetic electron fluxes are shown in Figure 9» The sharply peaked pro- 
files at 2, 5 and 10 keV agree very well with those given by Rees (1969). 

At 0.42 and 0.8 keV, however, the influence of elastic scattering becomes 
important and we find our profiles are broader and peak higher altitudes 
than predicted by Rees. 

A comparison of the ionization arising from electrons having ener- 
gies less than and greater than 500 eV is shown in Figure 10. Over a 
broad range of altitudes about 70 of the ionization comes from the inci- 
dent high energy electrons. Below 150 km, however, the ratio changes 
rapidly such that near and below the altitude of the peak ionization 
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rate, most of the ionization is caused by E < 500 eV electrons (including 
both degraded primaries and secondaries). 

The present ionization rates can be directly compared with, those 
computed from energy deposition curves and the assumption that 35 eV is 
expended for each electron-ion pair. The crosses in Figure 9 show the 
peak ionization rates derived from Fokker-Planck energy deposition cal- 
culations. These values show a systematic variation with incident energy, 
being too low at the lower energies and too high at the highest energy 
(10 keV). This behavior indicates the average energy needed to create 
electron- pairs is not constant with altitude but must depend to some ex- 
tent upon the energy of the incoming electrons. 

Finally, it is possible to give optical excitation rates for the 
monoenergetic fluxes. Figure 11 shows the computed 63OO A emission cor- 
rected for atmospheric quenching [k = 7*25 x 10~^n( Wgjsec - ' 1 "] . Scattering 
of the incident flux leads to a substantial upward flux which broadens 
the curves at high altitudes giving more emission than would be other- 
wise possible. Similar results are shown in Figure 12 for the 5577 A 
emission of atomic oxygen. 

Nighttime Aurora 

To illustrate a more practical situation, model computations have 
t-eeu ;«a.lr wr- t ng an incident electron flux spectrum measured during a 
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period of a moderate auroral activity at Fort Churchill, Canada. The 
results, described in further detail by Reasoner and Chappell (l973)> 
give a reasonably typical incident auroral spectrum having a peak flux 
in the vicinity of 11 keV. 

Using the incident flux, computations have been made of the elec- 
tron flux spectrum at various heights within the atmosphere. The results 
for the downward fluxes are shown in Figure 13 while the computed ioni- 
zation profile and several optical emissions are given in Figure lk. As 
was found for the monoenergetic cases, a pronounced minimum appears in 
the spectrums centered between -5 and 1 keV. Also, as the electrons 
penetrate more deeply into the atmosphere, the 11 keV peak is Attenuated 
and shifted towards lower energies. 

The influence of low energy incident particles can be seen in the 
ionization and optical emission profiles of Figure l4. Owing to differ- 
ences in the elastic scattering cross sections for 0 and N^, large changes 
occur in the low energy fluxes near 200 km, creating variations in the 
ionization rate and optical emissions. Below 200 km low energy electron 
transport becomes progressively more limited by collisions but above 
this height the electron free path soon becomes greater than the atmos- 
phere scale height so that not only are the incident low energy electrons 
reflected by atmospheric collisions but also there is an outward seepage 
of electrons from the lower altitude regions where substantial numbers 
of secondary electrons are produced. 



The ratio of dovmcoming to upgoing hemispherical fluxes is shown 
in Figure 15 as a function of energy. It appears that the albedo can 
range from 20$ to 100$. However, it is necessary to remember that the up- 
ward flux is composed of energy degraded and scattered primaries as well 
as secondary particles. 

Polar Cusp 

A more complete discussion of electron precipitation effects in the 
polar cusp will be given in a later paper. To illustrate the general 
effects we have adopted an incident power-law spectrum of the form 

cp = 2.34 x 109 cm -2 sec -1 ster -1 eV (l4) 

The downward flux spectrums resulting from this incident electron 
flux are shown for several altitudes in Figure 1 6 . As before, as one 
progresses downwards a minimum appears in the spectrums centered between 
250 to 700 eV. Similar results for E > 1 keV were obtained by Rees (l 969 ) 

The ionization profile and optical emissions corresponding to (l4) 
are shown in Figure 17. Because (l4) has been unrealistically extended 
to energies greater than 100 keV no peak is formed ^ in the ionization pro- 
file. The 6300 A emission does peak, however, as a consequence of quench- 
ing by N2. 
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Conclusions and Summary 


The present model of auroral electron behavior provides a substan- 
tially more comprehensive description than has been previously available. 

The computations at energies greater than 500 eV are similar in some ways 
to other models. However, the resulting description of the flux as a func- 
tion of energy, pitch angle and altitude is not used to compute an energy 
deposition rate, but rather is applied to ionization and excitation cross- 
sections to provide secondary and degraded primary production rates for 
electron energies in the range 0 - 500 eV. Using these production rates, 
the low energy program computes the remaining portions of the hemispherical 
flux spectrums in terms of upgoing and downgoing fluxes subject to various 
collisional processes. 

The results for the 2, 5 and 10 keV monoenergetic electrons agree 
well with the recent values obtained by Rees and Maeda (1973)* For 
these energies ionization deep in the atmosphere is not free to move sig- 
nificantly from its point of origin; i.e. although transport of secon- 
daries and degraded primaries is still present above 200 km, it cannot 
greatly alter the low altitude effects. 

At lower energies (l keV and below) the effects of particle trans- 
port become important and computations based upon energy deposition become 
inaccurate, especially when high altitude optical emissions are included. 

In such cases low energy particles may arise from both collisional degrada- 
tion and the presence of low energy electrons in the incident auroral 
spectrum. 
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Finally, the restriction of the present computations to hemispherical 
fluxes in the upwards and downwards directions could be removed using the 
methods of Stolarski (1972). However, until more detailed data are avail- 
able describing the angular factors associated with elastic and inelastic 
collisions, such an approach would probably not be worth the vastly in- 
creased computational effort. 
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FIGURE CAPTIONS 


1. Schematic diagram of the energy distribution of primary and secondary 
electrons created through impact ionization. 

2. Electron production rates as a function of energy in the range 0 to 500 eV 
at various altitudes for the nighttime auroral spectrum discussed in 

Section III. 

\ 

3* A comparison of discrete (solid curve) and continuous (dashed curve) 
energy loss rates for N^. The discrete losses include ionization and 
excitation while the continuous loss is based upon the Bethe formula 
of equation ( 3) • 

k. Electron production rates as a function of energy in the range 0 to 500 eV 
for a 10 keV Gaussian (monoenergetic) incident spectrum. The rates are 
normalized to unity incident (perpendicular) flux. 

5» Vertical profiles of downward and upward electron fluxes at selected 
energies for an incident Gaussian (monoenergetic) flux of 10 keV 
normalized to unity incident (perpendicular flux). 

6. Downward hemispherical fluxes at various altitudes for various Gaussian 

(monoenergetic) incident spectrums lying between 0.h2 and 10 keV. The 

% 

fluxes are normalized to unity incident (perpendicular) flux. 

7- Upward hemispherical fluxes at various altitudes for a 10 keV Gaussian 
(monoenergetic) incident spectrum. 
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8 . Ratio of §“/§ + for 10 keV Gaussian ( monoenergetic) incident electrons 
as a function of energy at several altitudes. Values less than 1 show 
where the upward flux is greater than the downward flux. 

9- Ionization rates per unit incident (perpendicular) flux for Gaussian 
(monoenergetic) fluxes in the range 0.42 to 10 keV. The symbol® 
indicates the results of energy deposition and the assumption that 
35 eV is expended for each electron- ion pair. 

10. Comparison of ionization produced by E < 500 eV electrons to the total 
ionization rate for different incident Gaussian (monoenergetic) fluxes. 

11. 63 OO A emission profiles (with quenching) for incident Gaussian (mono- 
energetic) spectrums normalized to unity incident (perpendicular) flux. 

12. 5577 A emission profiles for atonic oxygen with incident Gaussian (mono- 
energetic) spectrums normalized to unity incident (perpendicular) flux. 

13* Downward flux spectrums for a nighttime auroral spectrum obtained from 
rocket measurements by Reasoner and Chappell (1973)* For this example 
the hemispherical flux has been divided by 2tt ste radians. 

14. Optical emissions and the ionization profile associated with the incident 
nighttime auroral spectrum of Figure 13 . 

15. Ratio of §“/§ + as a function of energy at several altitudes for the night 
time auroral spectrum shown in Figure 13* 

16. Downward flux spectrums for a power law incident spectrum typical of the 
polar cusp. For this example the hemispherical flux has been divided by 
2rr ste radians. 

17. Opt. I cu i. emissions uml the ionlzutlon profile urising f rom the polur cusp 
type incident spectrum shown in Figure 16 . 
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TABLE 1 


MEAN ELECTRON EXCITATION PARAMETER, I 

l[eV}. 

n 2 8 o. 5 

0 2 92-0 

0 92.0 


TABLE 2 

SECONDARY ELECTRON ENERGY DEGRADATION COEFFICIENT 


N c 


0 , 


0 


E(eV) 

14.1 

18.9 

15-0 
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TABLE 3 


ADOPTED TOTAL ELECTRON IONIZATION CROSS SECTION PARAMETERS 



I 

A 

0 

V 

P 

V 

N 2 

19.96 

1.1 

O .85 

1 

1.2 

2 

°2 

18.5 

1.2 

0.87 

1 

1.1 

3 

0 

16.1 

0.84 

0.85 

0.3 

1.2 

l 
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TABLE 4 


ELECTRON IMPACT 

IONIZATION CROSS SECTION PARAMETERS 


State 

I 

A 

Q 

V 

y 

p 

X V 

15-58 

0.370 

0.80 

3.0 

1.0 

1.2 

aX 

16-73 

0.160 

O.83 

1.0 

1.0 

1.2 

B 2 £u + 

18.75 

0.073 

0.83 

2.0 

1.0 

1.2 

B2 "g 

22.0 

0.056 

O.83 

2.0 

1.0 

1.2 

cV 

23.6 

0.060 

O.83 

2.0 

1.0 

1.2 

Dissociative 

25.0 

O.38O 

O.96 

2.0 

1.0 

1.2 


12.1 

0.058 

0.80 

2.0 

1.0 

1.1 

a 4 fl u 

16.1 

0.150 

0.80 

2.0 

1.0 

1.1 

a2tt u 

16.9 

0.150 

0.80 

2.0 

1.0 

1.1 

bV. g - 

18.2 

0.130 

0.80 

2.0 

1.0 

1.1 

B 

23.0 

0.064 

0.80 

2.0 

1.0 

1.0 

Diss.( 0 +4 S) 

18.0 

0.400 

0.93 

3-0 

1.0 

1.1 

Diss.( 0 +2 D) 

22.0 

0.250 

0.93 

3-0 

1.0 

1.0 

4 s 

13-6 

0.290 

O.85 

1.0 

0.3 

■ 1.2 

2 D 

16.9 

O.36O 

0.85 

1.0 

0.3 

1.2 

2 P 

18.5 

0.190 

O.85 

1.0 

0.3 

1.2 
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TABLE 5 


ELECTRON IMPACT EXCITATION CROSS SECTION PARAMETERS 


State 

w 

(eV) 

A o 

Q 

V 

V 

A% + u 

6.14 

0.226 

3.0 

1.0 

1.0 

b \ 

7.30 

0.178 

3.0 

1.0 

3-o 

C \ 

11.03 

0.28 

3.0 

1.0 

3-0 

S 

9-10 

0.136 

1.0 

1.0 

1.0 

b \ 

12.85 

0.67 

0.75 

3 

1.0 

b V 

l4.o 

0-33 

0.75 

3 

1.0 

£ Rydberg 

13-75 

2.66 

0.75 

3 

1.0 

vib W = 10 

1.4 

1.5 x io5 

16 

9 

1.0 

g 

0.98 

0.0005 

3.0 

1.0 

3-o 

b V 

1.64 

0.0005 

3.0 

1.0 

3-0 

A& + 
u 

4-5 

0.021 

0.9 

1.0 

3-o 

B V 

8.4 

0.23 

0.75 

2.0 

1.0 

9-9 eV Allowed 

9-9 

0.68 

0.75 

3.0 

1.0 

E, Rydberg 

13-5 

2.77 

0.75 

3.0 

1.0 

vib 

-25 

9-57 X 10 -4 

.90 

1.0 

1.0 

1 D 

1.96 

0.01 

1.0 

2.0 

1.0 

1 S 

4.17 

0.0042 

1.0 

1.0 

0.5 

3s 

9-53 

0.0465 

1.0 

1.0 

0.5 

5 S 

9-15 

0.023 

2.0 

1.0 

1.0 

Z(hl=l, &s=0) 

14.2 

0.367 

0.75 

3.0 

1.0 

E(a s= 1) 

i4.7 

0.694 

2.0 

1.0 

1.0 

n(^l=0; &s=0) 

13-5 

0.043 

0.75 

1.0 

2.0 
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TABLE 6 


ADOPTED ELASTIC SCATTERING CROSS SECTIONS ( 10~ l6 cm 2 ) 


Electron 


Energy (eV) 

N 2 

°2 

0 

100 

8.1 

8.1 

4.9 

150 

7.2 

7.2 

4.4 

200 

6.3 

6.3 

3.8 

250 

5-5 

5-5 

3.3 

300 

4.7 

4.7 

2.9 

400 

3-5 

3-5 

2.1 

500 

2.4 

2.4 

1.4 
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